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Turbulent reacting flows in a generic swirl gas turbine combustor model are investigated both numerically and 


experimentally. In the investigation, an emphasis is placed upon the external flue gas recirculation, which is a 


promising technology for increasing the efficiency of the carbon capture and storage process, which, however, 


can change the combustion behaviour significantly. A further emphasis is placed upon the investigation of alter- 


native fuels such as biogas and syngas in comparison to the conventional natural gas. Flames are also investigated 


numerically using the open source CFD software OpenFOAM. In the numerical simulations, a laminar flamelet 


model based on mixture fraction and reaction progress variable is adopted. As turbulence model, the SST model is 


used within a URANS concept. Computational results are compared with the experimental data, where a fair 


agreement is observed. 
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Introduction 


In recent years, the link between the climate change 
and anthropogenic green-house emissions (especially 
carbon dioxide) has been well established. These pollu- 
tant emissions are mainly caused by the combustion of 
fossil fuels. One possibility for mitigating these pollutant 
emissions is the CCS technology [1]. Until now, the 
post-combustion capturing process based on amine 
scrubbing has been identified as a cost-effective and rea- 
lizable CCS technology [2]. The disadvantage of this 
technique is the significant deterioration of the overall 
plant efficiency, caused by the auxiliary power needed 
for the CO,-stripping process. This process can be im- 
proved by increasing the concentration of CO, in the 


Nomenclature 
C Reaction progress variable 
Cp Mean isobaric specific heat capacity - J/kgK 


k Turbulence kinetic energy — m’/s” 


exhaust gas, which can be achieved by external FGR [3]. 

The latter can change the combustion behavior of the 
gas turbine significantly. ElKady et al. [4], Winkler et al. 
[5] and Burdet et al. [6] reported an increase in CO emis- 
sions with increasing FGR ratios, which was essentially 
caused by a decrease in the combustion reactivity, due to 
the “dilution” of the flame by FGR, which, however, had, 
on the other hand the positive effect of reducing NO, 
emissions. The influence of hydrogen addition to NG was 
investigated by Winkler et al. [5]. Because of the in- 
creased reactivity, CO emissions at the equivalent FGR 
rates could be diminished, with only a slight increase in 
the NO, emission levels. 

As discussed in [1], a replacement of fossil fuels by 
renewable, alternative fuels is a further option for de 


b burnt 
st stoichiometric 
t turbulent 


lk Kolmogorov length scale - m 


p Static pressure - Pa 

P Probability density function 

Q Volumetric heat release rate — J/m’s 
t Time — s 

T Static temperature — K 

Xi Space coordinates — m 

Ui Velocity vector — m/s 

Y; Mass fraction of species i 

Z Mixture fraction 

Greek Symbols 

Bp Beta probability density function 

A Local finite volume cell size — m 

E Dissipation rate of k — m’/s° 

Oc Source term of the C transport equation 


Viscosity — kg/m's 
Kinematic viscosity — m?/s 


Flamelet parameter 


u 
v 
A 
P Density — kg/m? 
o Schmidt number 
9 


Variable to extract from flamelet libraries 
X Scalar dissipation rate — 1/s 
Subscripts 


ad adiabatic 


creasing CO, emissions. Kim et al. [7] and Lee et al. [8] 
investigated lean premixed combustion of LPG and SG 
blends in gas turbine combustors. Khalil et al. [9,10] in- 
vestigated the combustion characteristics of low calorific 
value gaseous and liquid fuels in premixed and diffusion 
mode. All these investigations were carried out using 
unvitiated combustion air. 

From the modelling point of view, there are two chal- 
lenges, namely, modelling of turbulence, and the model- 
ling of turbulent combustion, i.e. the interaction between 
combustion and turbulence. Turbulence modelling is an 
issue, due to the high degree of swirl that is normally 
encountered in gas turbine combustors, creating a 
strongly non-isotropic turbulence structure. Numerical 
investigations of isothermal swirling flows were, first, 
performed in [11-13]. However, these studies were per- 
formed for flows for relatively low Re (5000-7000) that 


do not represent typical gas turbine combustor conditions. 


Swirling flows with higher Re (Re >50,000) were inves- 
tigated in [14, 15]. The effect of the various turbulence 
modeling approaches based on URANS and LES proce- 


u unburnt 

Superscripts 

> Favre-averaged value 

> Reynolds-averaged value 


“ 


Favre fluctuational value 


Abbreviations 

BG Biogas 

CCS Carbon capture and storage 

EDC Eddy Dissipation Concept 

FGR Flue Gas Recirculation 

GI Grid Index 

GTC Gas Turbine Combustor/Combustion 
LFM Laminar Flamelet Method 

LES Large Eddy Simulation 

LPG Liquefied petroleum gas 

NG Natural Gas 

PDF Probability Density Function 

RANS Reynolds-Averaged Numerical Simulation 
RMS Root Mean Square 

SG Syngas 

SST Shear Stress Transport 

URANS Unsteady RANS 


dures was investigated in detail in [15]. In those studies, 
it was shown that the turbulence should be modelled 
within a three-dimensional and unsteady framework, 
which can also accommodate for the non-isotropic tur- 
bulence structure in large scales. These conditions are 
fulfilled best by the LES approach, and partially by the 
URANS approach. In the present analysis the URANS 
approach is employed. 

The further challenge for the simulation is the model- 
ling of turbulent combustion, i.e. the interaction of com- 
bustion reactions with flow turbulence. Many state-of- 
the-art turbulent combustion models are based on a glob- 
al reaction mechanism, assuming a small number of spe- 
cies, since, otherwise, the computational costs would 
strongly increase. However, in some cases, like the 
present one incorporating BG and SG fuels, a more de- 
tailed information on the chemistry and a consideration 
of a larger number of reacting species are required. A 
very efficient way of incorporating detailed chemistry 
without increasing computational costs is provided by the 
LFM. The use of the method for GTC have been rather 
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controversial, since it was claimed that the premixed 
combustion regime of GTC would not fall into the valid- 
ity region of the LFM. However, it was shown [16] that 
the LFM has a broader validity range than that, what was 
assumed based on purely dimensional grounds, and it can 
indeed lead to sufficiently accurate predictions [17] for 
GTC. In the present study, the LFM is adopted as the 
turbulent combustion model. 

Within the framework of the present investigation, the 
effects of external FGR on combustion behavior of natu- 
ral gas, typical biogases and syngas are investigated ex- 
perimentally and computationally. 


Experimental 


An atmospheric combustor was used, where the recir- 
culating flue gas was imitated by a natural gas fired con- 
densing boiler. By regulating the exhaust gas valve, a 
fraction of the boiler exhaust gas was diverted from the 
bypass conduit and mixed with fresh air. Upstream of the 
compressor, the flow rate as well as the oxidizer compo- 
sition was measured before being routed into the electric 
preheater, where the oxidizer is heated to typical com- 
bustor inlet temperatures. The partially premixed com- 
bustor is illustrated in Figure 1. 

The combustor consists of a radial swirler, a burner 
nozzle, an octagonal flame tube and an exhaust gas noz- 
zle. Four quartz glass windows are integrated into the 
flame tube for visual accessibility of the flame. The com- 
bustor is equipped with a pilot fuel system for ignition 
purposes as well as a main fuel system for main opera- 
tion mode. A detailed view of the injection system and 
the swirler is shown in Figure 2. 


windows exhaust gas 


Fig. 1 The model combustor. 
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Fig.2 Cross-sectional view of the combustor injection system 
(above) and the sketch of the radial swirler (below). 


Upon passing the combustor hood, the oxidizer stream 
enters the swirler radially and is redirected by the swirler 
vanes to impart a sufficiently high tangential velocity 
component to the flow. The swirler is made up of 12 
vanes with a constant vane angle of 45°, resulting in a 
swirl number of 0.93. Based on the local unburnt bulk 
axial velocity and the diameter at the exit of the con- 
verging-diverging burner nozzle, the Reynolds number 
turns out to be about 34,000. 

The combustor can be operated with NG and gas mix- 
tures. The gas mixing device made up of eight separate 
mass flow controllers can continuously provide a mixture 
of On, N2, CO, CO, CHy, Hz and C3Hg in a required 
composition. For NG operation - the pilot fuel system is 
always operated by NG - similar mass flow controllers 
are connected to the NG grid. The exhaust gas tempera- 
ture is measured with a Type S thermocouple that is in- 
corporated into a ceramic exhaust gas probe. The compo- 
sition of the oxidizer and emissions are measured by an 
industrial system including sensors for O2, CO2, CO, NO, 
and unburnt hydrocarbons (UHC). The measurement 
error was less than or equal to 2% of the measured value 
for CO, NO, and UHC analyzers, and less than or equal 
to 1% of the full scale value for the O, and CO, analyzers. 
The detection limits for CO and NO, emissions were 
specified as 0.4 ppm and 0.3 ppm, respectively. A more 
detailed description of the rig is given in [18]. 

The oxidizer compositions that are considered in the 
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present study are shown in Table 1. A list of the investi- 
gated gaseous fuels is provided in Table 2. In the present 
work, NG, BG5 and SG fuels are considered. For all ex- 
periments, the combustor inlet temperature is kept at the 
constant value of 573 K. The adiabatic flame tempera- 
tures are kept in the ranges 1500 K-1550 K, 1600 K - 
1650 K and 1450 K-1500 K for the NG, BG, and SG fu- 
els, respectively. Lower adiabatic flame temperatures 
were chosen for the SG fuel for preventing the risk of 
flashback [19], due to its high hydrogen content (Table 2). 


Table 1 Composition of investigated oxidizers. 


FGR oxidizer fraction in vol % 
% O, CO: H0 Ar N2 
0 20.6 0 1 0.9 715 
20 17.9 2.2 2.3 0.9 76.7 
40 13.5 5.9 3.6 0.9 76.1 


Table 2 Composition of investigated gaseous fuels. 


fuel fraction in vol % 


Gas 
CH; CO, G&H CO H O No 
BGI 65 30 : : : 2 3 
BG2 55 40 - z : 2 3 
BG3 45 50 - : = 2 3 
BG4 58.5 27 - - 10 18 27 
BGS 495 36 - : 10 18 27 
BG6 405 45 4 : 10 18 27 
SG 10 22 4 22 40 - 2 
NG 925 13 5.2 : : - 1 
Modelling 


A block structured grid, comprising 1.2 million finite 
volumes is used. Views of the surface grid are illustrated 
in Figure 3. At the inlets, constant profiles are prescribed 
for velocities, temperature and mass fractions that result 
from the combustor data. At the outlet, fully-developed 
boundary conditions are applied for convective-diffusi- 
vely transported variables, along with a radial equili- 
brium for the static pressure. Although there exist more 
sophisticated formulations such as nonreflecting boun- 
dary conditions [20] they are not implemented, as the 
above-mentioned formulation is found to be sufficiently 
adequate for the present low Mach number case. The 
no-slip walls are assumed to be adiabatic. 

For the computational analysis, the open source finite 
volume CFD code OpenFOAM [21] is used. The pres- 
sure-correction scheme PIMPLE is applied for the veloc- 
ity-pressure coupling [21]. As turbulence model, the 
Shear Stress Transport (SST) model [21] is used within a 
URANS framework. A gradient-diffusion approximation 
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is used for scalars, assuming constant turbulent Schmidt 
numbers of value 0.7. 


Fig.3 Views ofthe surface grid. 


A second order upwind discretization scheme is ap- 
plied to the convective terms for all variables. A first or- 
der Euler scheme is used for time discretization, as sta- 
bility problems were encountered using a second order 
scheme. The time step size is chosen to ensure cell Cou- 
rant numbers to be smaller than unity. Starting from an 
initial field, simulations were performed for a period long 
enough to allow the establishment of a quasi-periodic 
flow field that is no more depending on the initial condi- 
tions. After this state, the time-averaging was started, and 
continued until the time-averaged fields did not show any 
remarkable change in time. 

In the present analysis, the Laminar Flamelet Method 
(LFM) is applied as turbulent combustion model. Ac- 
cording to the assumptions of the LFM [22], for a steady, 
1D, adiabatic, laminar diffusion flame, all thermo- 
chemical dependent variables (g) can be expressed as 
functions of the mixture fraction Z and the stoichiometric 
value of the scalar dissipation rate 7,,as 

9=0AZ,%4) (1) 

The mixture fraction Z denotes the mass fraction of 
fuel stream locally in the unburnt mixture. For Xs, quite 
commonly, an assumed functional dependence on Z (e.g. 
as the one suggested by Peters [22]) is used. Such func- 
tional relationships obtained by performing 1D laminar 
flame calculations based on 


2 
xOF 
Pa oz =m; =0 (2a) 
x= 2D,(VzZy (2b) 


once for all, and made available in tabulated form 
(flamelet libraries) for the subsequent CFD analysis (In 
Eq. 2, m; and D, denote the source term and the diffu- 


sion coefficient, respectively). Assuming Z and Xs are 


independent, based on presumed PDFs that are to be ob- 
tained from the results of the turbulent flow, the average 
thermo-chemical dependent variables can, then, be ob- 
tained from (the expression is written for a Fa- 
vre-averaged quantity, but holds analogously for a Rey- 
nolds-averaged quantity) 


co pl 
B= | | eZ. x)PZ)Pleu)dZd zx, 63) 
Typically, a single-delta PDF is assumed for y,. The 
presumed PDF for Z is typically controlled by its two 
moments, i.e. the average value Z and the variance 


Z"? . The modelled time-averaged differential transport 
equations of these variables [23] are: 


az A A 
op aZ) offe uZ] ua 
ot Ox; Ox; |o oO, J Ox; 


opZ”  pii,2"”) ô (eZ 
ôt ôx; Ox ; j i Ox; (4b) 
-2% OZ OZ a7" 
o, Ox; Ox; k 


c being a model constant (c=2.0). The above outlined 
LFM is known to perform well for a large class of 
non-premixed flames. However, it is known [22] that it 
cannot accurately describe phenomena such as local ex- 
tinction, re-ignition, flame lift-off. Thus, it cannot be 
assumed to be convenient for premixed/partially pre- 
mixed flames such as the currently considered one. 
Therefore, in the present work an alternative flamelet 
model is adopted, which can better cope with premixed/ 
partially premixed combustion. The model was proposed 
by Pierce and Moin [24] within an LES framework, who 
parametrized the flamelets based on a new flamelet vari- 
able 4 that is defined through the so-called reaction 
progress variable C as given by 
A = max(C) (5a) 
The reaction progress variable can be defined in dif- 
ferent ways. In the present work, a temperature-based 
definition is preferred (the original species-based defini- 
tion is not convenient in the present case with FGR) 
= T-T, u 
T,-T, 


u 


(5b) 


Hence, the laminar flamelet functional relationships 
are established as 
y= 9(Z.a) (6) 
Thus, the average values thermo-chemical variables 
are obtained in a similar fashion to Eq. (3). In the present 
work, a single-delta PDF is assumed for 4. For Z, a beta 
PDF (P) is used. Thus, the average values are obtained 
from 


Ge [, f o(z,a)alz;Z,2" \a(a-Z)azaa (7) 


As å is connected to C, for obtaining the correspond- 
ing information, a differential transport equation is 
solved for C, derived following the Bray’s approximation 
[25] of the chemical source term reading as 

ope olp č) _ a [s n) a + 
+ + =0c (8) 
Ot Ox; Ox; |\O o, J Ox; 
In Eq. (8) Oc is defined as 


bc = — g 
* eT, —T,) ý 

where Q and cp denote volumetric heat release rate 
and the mean isobaric heat capacity, respectively. The 
volumetric heat release rate Q is obtained from the 1D 
laminar flame libraries. 

The flamelet libraries are constructed by 1D steady, 
adiabatic laminar flame calculations using the Flame- 
Master code [26], a priori to the field calculations of the 
turbulent reacting flow by means of CFD. As the under- 
lying reaction mechanism, the GRI Mech 3.0 [27] is used, 
assuming a Lewis number of unity for all species. 


Results 


Assessment of the Computational Grid 


In the computational study, a formal grid indepen- 
dence study was not performed. An analysis on the tur- 
bulent scales implies, however, that the present grid res- 
olution is sufficiently fine. Within the framework of LES, 
different measures were proposed for assessing the grid 
resolution, such as the Grid Index (GI) defined as ratio of 
the local grid size (assumed to be given by the third root 
of the cell volume) to the Kolmogorov length scale. 


cia (10) 
lk 


1 
ad 
y= Bi a1) 


According to Celik et al. [28], GI should be smaller 
than 25, for achieving sufficient accuracy for LES. In the 
present work, not an LES, but a URANS approach is 
used. Still, the LES grid resolution criteria can be re- 
garded to be useful to have an indication of the grid res- 
olution quality. Figure 4 presents the distribution of GI 
calculated for the NG flame with 20% FGR, at a time 
step, in a plane through the swirler (at the channel mid- 
height), which is a critical region for the mixing and where 
the length scales are small. One can see that the values do 
not exceed 25 considerably and remain lower than 50. 
Thus, considering that not an LES but a URANS ap- 
proach is used in the present study, which is theoretically 
less demanding as far as the grid resolution is concerned, 
it can be assumed that the present grid adequately fine. 


with 
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Fig. 4 GI distribution at an instant of time, in a plane through 
the swirler (at the channel mid-height) for the NG 
flame with 20% FGR. 


Natural Gas Flames 


Calculated distribution of instantaneous velocity mag- 
nitude, time-averaged streamlines, iso-surface of high 
vorticity and distribution of instantaneous and time- av- 
eraged temperature are displayed in Figure 5. It can be 
observed that a large inner recirculation zone is formed at 
the exit of the burner nozzle, at the frontal stagnation of 
which, the flame is anchored. The unsteady features of 
the predicted flow field exhibiting a precessing vortex 
core can also be observed. 

As it can be seen in the predictions presented in Figure 
5, the flame is positioned at the exit of the burner nozzle, 
as confirmed by the experiments. The combustion takes 
place in the premixed mode, where maximum tempera- 
ture of about 1500 K are achieved behind the slightly 
curved flame brush (Fig. 5e). 

The measured and predicted temperature variations 
along the combustor length are presented in Figure 6 for 
the NG flame with 20% FGR, along the combustor 
length (axial distance is measured from the position of 
the burner nozzle outlet (Fig. 1)), and along a line paral- 
lel to axis at a radial position r = 30 mm (This line 
touches the burner nozzle outlet edge (Fig. 1)). Both ex- 
periments and predictions show a rapid temperature rise 
across the rather thin flame front, as already indicated by 
the temperature plot of Figure 5. Since the flame front is 
curved and the flame brush is comparably thicker near 
the burner nozzle edge (Fig. 5) the temperature increase 
observed along the r = 30 mm line (Fig. 6b) is less steep 
compared to the temperature rise along the combustor 
axis (Fig. 6a). The predicted rise of temperature in the 
flame front shows a good agreement with the experi- 
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ments (Fig. 6). The burnt gas temperature is overpre- 
dicted by about 100-150 K (Fig. 6). This may be due to 
the assumption of an adiabatic flame in the predictions. 
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Fig. 5 Predictions for the NG flame with 20% FGR (a) Veloc- 
ity magnitude in combustor mid-plane at an instant of 
time, (b) time-averaged streamlines in combustor mid- 
plane, (c) iso-surface of high vorticity at a time step, (d) 
temperature field in combustor mid-plane at an instant 
of time, (e) time-averaged temperature field in com- 
bustor mid-plane. 


Measured and predicted mass fractions (in ppm @15% 
O2) of CO and NO at the combustor exit are presented in 
Figure 7, for the NG flame with 0% FGR and 40 % FGR. 
The predicted emission mass fractions agree with the 
measured values especially good for the 40% FGR (Fig. 7). 

For 0% FGR, the agreement for the NO mass fraction 
is better compared to CO mass fractions, whereas for the 
latter, the calculations clearly overpredict the very low 
measured value. The overall agreement can still be con- 
sidered to be quite satisfactory (Fig. 7). An interesting, 
and, having the thermal NO formation mechanism [29] in 
mind, rather unexpected trend is that the NO emission 
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Fig. 6 Time-averaged measured and predicted temperature 
variation along combustor length for the NG flame 
with 20% FGR, (a) along burner centerline, (b) along 
r= 30 mm line. 
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Fig. 7 Measured and predicted mass fractions (in ppm @ 
15 % Oy) at the combustor exit for the NG flame with 
0% FGR and 40% FGR, (a) CO mass fractions (b) NO 
mass fractions. 


increases with increasing FGR. One would normally ex- 
pect the contrary since both major contributors to NO 
formation rate, i.e. the oxygen concentration and temper- 
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ature (thermal NO), decrease with increasing FGR. This 
trend is clearly observed in the measurements and only 
slightly predicted by the calculations. This trend will be 
discussed in more detail in the subsequent section on BG 
flames. 


Biogas Flames 


Predicted time-averaged temperature fields for the BG 
flame (the BGS fuel) for different levels of FGR (0% 
FGR, 20% FGR, 40% FGR) are presented in Figure 8 in 
combustor mid-plane. It is interesting to note that flue 
gas recirculation has a very remarkable influence on the 
thickness of the flame brush according to the predictions, 
as observed especially for the 40% FGR case. Further- 
more, it can be observed (Fig. 8) that the BG flame exhi- 
bits a generically more strongly curved/inclined flame 
front compared to the NG flame, which exhibits a com- 
parably flat one (Fig. 5e). The higher hydrogen content of 
the BG fuel may be playing a role here that allows a 
higher penetration of the flame in the core region, due to 
a higher flame speed, whereas the boundary regions may 
be protected against flame propagation due to the locally 
higher strain rates. 

Measured and predicted mass fractions (in ppm @15% 
O2) of CO and NO at the combustor exit are presented in 
Figure 9, for the BG flame with 0% FGR, 20 % FGR and 
40 % FGR. One can observe that both the measured CO 
and NO emissions for the BG flame steadily increase 
with increasing FGR ratio, which is predicted qualita- 
tively well by the calculations (Fig. 9). As mentioned 
above, it is rather unexpected that NO emissions increase 
with increasing FGR ratio, if one considers the thermal 
NO formation path. The reason for the observed behavior 
is that increasing FGR ratios beyond a certain level (10% 
-15% FGR) shifts the operation point nearer to stoichi- 
ometry, which is beneficial for another NO formation 
path, namely for the prompt NO via the Fenimore me- 
chanism [29], which is dependent on hydrocarbon radi- 
cals and occurs mainly during the heat release. The latter 
is prolonged with increasing FGR ratios as manifested by 
the predicted thickening of the flame brush with increas- 
ing FGR ratio (Fig. 8). Similar trends, i.e. increasing NO 
emissions with increasing FGR ratio were also observed 
in the previous studies of other authors [30]. At this stage, 
it should be pointed out that these effects can only be 
predicted by a model incorporating detailed chemistry, 
and the presently employed LFM represents a very effi- 
cient alternative in that respect. On the other hand, it 
shall also be pointed out that the vertical scale of the CO 
plot (Fig. 9a) is a logarithmic one. Thus, the quantitative 
disagreement between the measurements and predictions 
is actually larger than it seems to be at the first sight. The 
comparably inferior quantitative prediction for CO may 
be due to fast chemistry assumption underlying the LFM 
being not very suitable for representing the CO kinetics. 
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Fig. 8 Predicted time-averaged temperature fields in combus- 
tor mid-plane for the BG flame, (a) 0% FGR, (b) 20% 
FGR, (c) 40% FGR. 
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Fig. 9 Measured and predicted mass fractions (in ppm @ 
15 % O,) at the combustor exit for the BG flame with 
0% FGR, 20% FGR, and 40% FGR (a) CO mass frac- 
tions (b) NO mass fractions. 


Syngas 

Predicted time-averaged temperature fields in com- 
bustor mid-plane are presented for the SG flame for 20% 
FGR and 40% FGR in Figure 10. One can observe (Fig. 
10) that the inclination of the flame front is even stronger 
for the SG flame, compared to the BG flame (Fig. 8b). 
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This may, again, be attributed to the even higher hydro- 
gen content of the SG (Table 2). The thickening of the 
flame brush (Fig. 10) with increasing FGR ratio is quite 
moderate compared to the BG flame (Fig. 8). 

Measured and predicted mass fractions (in ppm @15% 
O2) of CO and NO at the combustor exit are presented in 
Figure 11, for the SG flame for 20 % FGR and 40 % 
FGR. It can again be observed that the calculations pro- 
vide a qualitatively well prediction of the experiments. 
The quantitative agreement is better and rather pretty 
good for the NO mass fractions (Fig. 11b), whereas the 
quantitative discrepancies are larger for the CO mass 
fractions (Fig. 11a). 
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Fig. 10 Predicted time-averaged temperature fields in com- 
bustor mid-plane for the SG flame, (a) 20% FGR, (b) 
40% FGR. 
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Fig. 11 Measured and predicted mass fractions (in ppm @ 15% 
O2) at the combustor exit for the SG flame with 20% 
FGR and 40% FGR (a) CO mass fractions (b) NO 
mass fractions. 


Conclusions 


Turbulent flames in a generic swirl gas turbine combustor 
are investigated numerically. The turbulent flow is mod- 
elled within a URANS framework, using the SST turbu- 
lence model. For modeling the turbulence-chemistry in- 
teraction, a laminar flamelet model based on the mixture 
fraction and reaction progress variable is used, coupled 
with a presumed probability density function approach. 
Natural gas, biogas and syngas flames with and without 
external flue gas recirculation are investigated. The nu- 
merical results are observed to show a fair agreement 
with the experimental data. The slight overprediction of 
the combustor exit temperature is assumed to be caused 
by the assumption of adiabatic flame in the mathematical 
model. The model will be improved in the future to in- 
clude non-adiabatic effects, which is also expected to 
lead to a more accurate prediction of NO emissions. Ac- 
curacy of the CO emissions were in general not as good 
as that of NO emissions. This may be due to the applied 
laminar flamelet model, which may not be very suitable 
for the relatively slow CO chemistry. Alternative turbu- 
lent combustion models will also be explored in the fu- 
ture. 
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